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y—i '■ ABSTRACT 

' Based on the CCD observations of remote young open clusters NGC 2383, NGC 2384, 

t> ! NGC 4103, NGC 4755, NGC 7510 and Hogg 15, we construct their observed lumi- 

nosity functions (LFs). The observed LFs are corrected for field star contamination 
determined with the help of galactic star count model. In the case of Hogg 15 and 
NGC 2383 we also considered the additional contamination from neighboring clusters 
NGC 4609 and NGC 2384 respectively. These corrections provided the realistic pat- 
| tern of cluster LF in the vicinity of the MS turn on point and at fainter magnitudes, 

"--»« i revealed the so called H- feature arising due to transition of the Pre-MS phase to MS, 

*-£h ' which is dependent on the cluster age. The theoretical LFs were constructed represent- 

ing a cluster population model with continuous star formation for a short time scale 
O ' an d a power law Initial Mass Function (IMF) and these were fitted to the observed LF. 

| As a result we are able to determine for each cluster a set of parameters, describing 

c/3 , cluster population (the age, duration of star formation, IMF slope, and percentage of 

field star contamination). It was found that in spite of the non-monotonic behavior of 
observed LFs, cluster IMFs can be described as the power law functions with slopes 
• »h . similar to Salpeter's value. The present Main Sequence turn on cluster ages are several 

y\ ' times lower than those derived from the fitting of theoretical isochrones to turn off 

^ , region of the upper Main Sequences. 

Key words: keywordl: specificationl - keyword2: specification - ... (see 4th issue 
of each volume) 
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1 INTRODUCTION 

Young clusters are natural laboratories for the study of vari- 
ous issues related to star formation processes, e.g., the initial 
mass function (IMF), duration of star formation etc. In 
recent years, some remote young clusters have been studied 
using BVRI CCD obse rv ations ( e.g., see lSagar fc Griffithsl 
19911. IPhefps fc Janes] Il993t ISaear fc Cannon! Il995l 
Sagar fc CannorJll997l ISaear fc Griffithslll99ct I >aear et all 

Il999l. ISanner et alJ l200lL 
Pandev et alJl200ll iPiatti et alJ 120021 . These observations 
were used for the construction of color-color diagrams, 
color-magnitude diagrams (CMDs) and luminosity func- 
tions (LFs). These diagrams were used for the determination 
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of cluster reddening, distance and age while the LFs were 
transformed to stellar mass spectrum, which can be identi- 
fied with the IMF for young (age ^ 100 Myr) star clusters. 
In these conventional studies of young star clusters, cluster 
age, LFs and MFs could not be determined accurately due 
to reasons discussed below. 

Firstly, for reliable separation of cluster members from 
the foreground and background field stars present in the di- 
rection of distant (d > 1 kpc) clusters, accurate kinematic 
(proper motion and/or radial velocity) measurements are es- 
sential. Unfortunately, even present day most accurate Hip- 
parcos proper motion data (available for brightest V < 13 
mag stars only) are not sufficient for this purpose. In such 
circumstances the field star contamination is removed statis- 
tically by observing comparison areas (blank fields) adjacent 
to the cluster field. It should be noted, however, that in most 
cases of the conventional general-purpose CCD photometry, 
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which are not specially intended to e.g. the IMF study, such 
regions are not observed, and hence even statistical means 
of the data refining could not be used. This leads to a loss 
of information suitable otherwise for the LF construction. 

The other problem in the CMDs of young star clus- 
ter is the precise determination of the Main Sequence (MS) 
turn off point, which is used for age determination, since 
for young clusters due to their small age, the upper MS is 
practically unevolved. Consequently, it becomes difficult to 
evaluate the accurate cluster age, and instead an upper esti- 
mate of the cluster age is provided. The difficulty is further 
increased by effects like stellar rotation and binarity, and by 
the morphology of the upper MS, which is practically verti- 
cal for early type stars in the colours (B — V), (V — R), or 
(V-I). 

These problems can be avoided if the cluster age is esti- 
mated using the location of the MS turn on point, the place 
where Pre-MS stars join the MS in the cluster CMD. For a 
young coeval cluster this point resides about 7 mag fainter 
the brightest termination point of the MS. As a result of 
the deep CCD observations of the clusters, this region can 
be seen in the CMD, and hence is available for the anal- 
ysis. Many of the above effects (except that of unresolved 
binaries) are much weaker in the turn-on point region, which 
makes the turn-on dating method to be more attractive than 
the turn off one. There are, however, practical problems in 
locating the turn-on point e.g., presence of strong field star 
contamination hides it more strongly than the brightest MS. 
The other limitation of the turn-on method is cluster age. 
As shown by evolutionary calculations, only for extremely 
young clusters with ages of order of a few Myr, the Pre-MS 
branch is raised sufficiently above the ZAMS and could be 
easily recognized in the CMD. For older clusters the Pre-MS 
branch deviation from the ZAMS diminishes with increasing 
age and could be seen better as a detail of the LF (as it is 
discussed in Sec. At cluster ages exceeding that of the 
Pleiades, the Pre-MS branch could not be identified with 
confidence both in the CM D and in the LF and the turn on 
method cannot be applied feelikov et al.lH998l) . 

The other related issue is the fine structure of the stel- 
lar LF, located in the vicinity of the MS turn on point, 
which can mask the IMF shape if one does not take it 
into account. This feature can produce a false flattening or 
even depletion in the LF, which are frequently considered as 
an evidence of the flatten ing /turning over of the IMF (see 
IPiskunov fc Belikovlll996l for references) . This in turn might 
have an important consequence with respect to the IMF uni- 
versality and other similar issues. Since luminosity function 
of the turn on point depends on cluster age, the position of 
the LF detail also varies with time, and may be used as a 
kind of a standard candle for age de termination from the LF 
analysis feelikov fc Piskunovlll997h . 

Thus, the MS turn on region of young cluster CMDs is 
very important for cluster dating, or should be taken into 
account when one is analyzing the LF. It can be easily iden- 
tified, however, only in the case of a few selected clusters. 
The well known examples of such clusters are NGC 2264, 
NGC 6530, ONC, and some others. Normally young and 
especially remote clusters, which are buried in rich fore- 
ground/background show, neither distinct turn ons nor PMS 
branches. 

The aim of present study is to reveal the MS turn on 



point information hidden in the existing CCD observational 
data of 6 young open cluster using a new approach, which 
does not require observations of blank fields. It should be 
noted, that these clusters were subject of a conventional 
study, and many of them already have age and IMF deter- 
mination. Unlike these studies, we will not convert observed 
LF to the IMF, but in order to avoid the above mentioned 
problems, we will construct theoretical LF, which should 
properly reflect the behavior of the observed ones in the 
vicinity of the MS turn on point. Wherever the field star 
contaminations could not be determined using the observa- 
tions of an offset field region, they are estimated using galac- 
tic model for star counts for the surroundings of the cluster 
under study. Then the theoretical LF will be fitted to the 
observations by varying the star formation parameters de- 
scribing actual stellar population (cluster age, star formation 
duration, IMF slope, a percentage of observed field stars). 
Considering its success, we hope that this approach can also 
be implemented to the other young open clusters having 
similar or more deep CCD data, w hich become presently 
available (see e.g. iKalirai et aHl200lb . 

In Section [5] we detail theoretical approach to LF of a 
young cluster and consider its fine structure in the context 
of current study. In the Sec. |21 we summarize the data used 
in the present work. Sec. [I] describes construction of the 
observed and theoretical LFs. We consider here the major 
effects which should be taken into consideration, define the 
model of cluster population and its parameters, and describe 
the fitting procedure. In Sec.|^|we discuss the derived results, 
and summarize them in Sec. |H] 



2 THEORETICAL APPROACH 

Since stellar mass m can be directly measured in very rare 
and specific cases, the IMF is not observed directly, but is 
converted from the observed distribution of stars over their 
absolute magnitudes dN/dM, called LF 4>{M). For a cluster 
of age t (in years throughout the paper if not mentioned 
otherwise), the LF <f>t[M) is related to the IMF f(m) via 
time-dependent mass-luminosity relation (MLR) m t (M) as 
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Since cluster stars evolve off or approach to the MS with 
a rate which is dependent on their mass, apart from their 
stay on the MS, the cluster MLR should also evolve with 
time, being different at every moment from the standard 
one (usually adopted to be the MLR of the MS stars) . This 
deviation is especially strong for red giant and Pre-MS stars. 
Besides, the "instant" MLR has a definite fine structure due 
to presence of quasi-horizontal Post-MS and Pre-MS por- 
tions of corresponding isochrone. The structure itself is not 
too prominent to influence the LF, but it produces bumps 
and dips in the MLR derivative, and hence in the LF. The 
strongest fluctuation of the derivative occurs near the MS 
turn off point, but due to the small number of stars observed 
there they could not be easily distinguished from statistical 
fluctuations of the LF. In contrast, in the vicinity of MS turn 
on point, a sufficient number of stars coupled with strong 
MLR derivative bumping produce a feature, which can be 
observed in LFs of star clusters, which are young enough 
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Figure 1. Theoretical color-magnitude diagram (Panel a) and luminosity function (Panel b) of a young open cluster. Solid curves are for 
an isochrone of log t = 7.0, including both Pre-MS and Post-MS evolutionary stages, and corresponding luminosity function, constructed 
from the isochrone. T he dotted curves are th e ZAMS and "Main Sequence" luminosity function. Thin curves in panel (a) are Pre-MS 
evolutionary tracks of iPalla fc Stahlerl ll993fl . The crosses and adjacent figures separate different evolutionary stages as it is explained 
in the text. The enlarged MS turn on point area is shown in the panel encapsulated in the CMD. 



to display a branch of pre-M S stars. This detail was called 
bv lPiskunov fc Belikovl 1^996) as H-feature, as it appears at 
the location, where the hydrogen burning starts in the cores 
of the pre-MS stars. 

In order to get an idea on general structure of the LFs of 
young star clusters and their relation to main evolutionary 
stages of stars, we show in FigQthe CMD together with an 
isochrone of logt = 7.0, and corresponding LF constructed 
according to eq. Q with help of Salpeter IMF (f{m) oc 
to -1 ' 35 ). The isochrone was constructed from a set of stellar 
models described in Sec. 14.41 Our calculations show, that in 
general this behavior of the LF does not differ for various 
sets of models, but details (the size and amplitude of the LF 
fragments) may be somewhat different. 

As one can see from Fig. 0for young clusters the LF is 
non monotonic even for monotonic IMF, and consists of a 
number of monotonic portions related to different evolution- 
ary stages. Comparison with evolutionary tracks shows that 
segment (0,1) of the LF corresponds to convective portion of 
a Pre-MS track, segment (1,2) to radiative one, (2,3) can be 
identified with MS turn on point in the CMD, and segment 
(3,4) belongs to the Main Sequence. Due to non-monotonic 
behavior of the isochrone near MS turn on point, the seg- 
ment (2,3) is in fact a superposition of MS and Pre-MS 
stellar evolutionary stages and pure MS starts at somewhat 
brighter magnitudes. The last segment (4,5) corresponds to 



Post-MS stages of stellar evolution. We call hereafter for 
convenience purposes segment (1,2) as the LF H-dip, and 
segment (2,3) as LF H-maximum. F or comparison in F ig. U 
we show position of ZAMS models of lTout et al.l ( 1996), and 
ZAMS-LF constructed from these models. Good agreement 
between ZAMS- and LFs is observed between both turn-over 
points, but beyond MS both LFs differ considerably. 

Since the absolute magnitude of a cluster MS turn on 
point evolves with time, the H-feature should move faint 
ward with increasing of cluster age. The feature not only 
moves to fainter magnitudes, but also degrad es with time. 
As calculations of iBelikov fc Piskunovl il997f) show it dis- 
appears completely for logt ^ 8.2. This puts a limit on the 
H-calibration applicability. Again the time scale of H-feature 
existence depends on specific models, and can be used for 
Pre-MS model verification. In principle this makes the LF 
of young clusters a good tool for cluster age determination. 

An immediate conclusion, which could be drawn from 
the above analysis is that the existence of any bumps/gaps 
in the LF at or near the MS turn on point of young clusters 
should not be regarded as an evidence of non-monotonous 
behavior of the IMF. The existence of the Pre-MS detail 
in the LF may also result in the mis-interpretation of mass 
spectra of young open clusters. When the limiting magni- 
tude of a survey falls within the radiative dip one could ob- 
serve the LF turn-over, which can be wrongly interpreted as 
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Table 1. Specification of studied clusters and used observation data 



Cluster 


IAU 


d' 


Area 


Number 


Vlim 




V — My 


E{B - V) 


log* 


Refe- 




number 




(□') 


of stars 


(mag) 


(mag) 


(mag) 


(mag) 




rence 


NGC 2383 


C0722-208 


5 


31.8 


588 


21.8 


17.6 


13.3 


0.22 


8.45-8.6 


[1] 


NGC 2384 


C0722-209 


5 


27.8 


256 


20.9 


18.9 


13.2 


0.22-0.28 


7.1-7.3 


[1] 


NGC 4103 


C1204-609 


5 


19.4 


176 


19.8 


17.3 


12.5 


0.31 


7.5 


[2] 


Hogg 15 


C1240-628 


2 


6.9 


337 


21.7 


19.6 


16.5 


1.15* 


6.78 


[3] 


NGC 4755 


CI 250-600 


10 


38.8 


576 


20.2 


18.1 


12.9 


0.41* 


7.0 


[■1] 


NGC 7510 


C2309+603 


6 


38.4 


423 


21.1 


16.8 


16.0 


1.12* 


<7.0 


[5] 
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* non uniform extinction 



a consequence of data incompleteness or the IMF turnover. 
Similarly the steep slope of convective portion of the LF is 
not an evidence of the IMF steepening. 



3 OBSERVATIONAL DATA 

We have used CCD photometric data of young (age sev- 
eral tens of Myr) remote and compact clusters presented in 
papers listed in Table Q (hereafter referred as original pa- 
pers). The small areas occupied by these clusters and their 
sufficiently deep photometry provide high degree of data 
completeness, which is very important for LF construction. 
However, their large distance coupled with lack of kinematic 
data prevent the separation of cluster members from the 
field stars and requires special approach to tackle this issue. 

In Table Q the parameters of the clusters under study 
are listed. Cluster designations are shown in columns 1 and 
2, in colum n 3 we list v alues of cluster angular diameters 
taken from iLvngal <ll987l) . The most of the other parame- 
ters are taken from the original papers mentioned in the 
last column. Columns 4 through 9 contain area covered by 
CCD frames, number of stars within the frames, limiting 
and completeness ^/-magnitudes, apparent distance moduli 
and reddening values. Color excesses, marked with asterisks 
represent average values in variable extinction fields. In col- 
umn 10 we show log(age) as it was determined in referred 
studies. We define limiting magnitudes of a survey as faintest 
magnitude of a star in cluster sample. Data completeness in 
magnitude are estimated following the procedure described 
in Sec. IQ1 

We compared original data of TableHlwith lLvngal dl987l) 
catalogue and recent list of reddening values, distance mod- 
uli and ages, provided by Loktin 1 (private communication, 
referred hereafter as LGM2.2). We find that the extinction 
values from Tableware in excellent agreement with the data 
of both catalogues and the average color excesses from differ- 
ent lists agree within a few hundredth of magnitude, except 
in the case of NGC 7510. In this case, the extinction data as 
in LGM2.2 differ from that of TableQby 0.26 mag. Original 
cluster ages also show reasonable agreement (with average 
spread of the order of a few tenths in logt) for all c l usters 
except NGC 2383, which according to both ILvngal dl987T) 
and LGM2.2 has an age of a few tens of Myr. 

1 Loktin A. V., Gerasimenko T.P., Malysheva L.K., 2001, Homo- 
geneous Catalogue of Open Cluster Parameters, Release 2.2. See 
also WEBDA database at URL http:/ /obswww.unige.ch/webda/ 



Unfortunately such a good agreement is n ot observed in 
the case of cluster distances. Comparison with lLvngal (I1987T) 
referred distance estimations shows good agreement for all 
clusters, except NGC 2383 and NGC 2384, where data of 
TableQare overestimated by about 1 mag. The new scale of 
LGM2.2 indicates even stronger (up to 1.5 mag) disagree- 
ment for these and two more (Hogg 15, and NGC 7510) 
clusters. Though these disagreements are unpleasant, it is 
hardly avoidable as photometric distances of remote clusters 
are derived from a very steep upper MS. In present study we 
will use original distances, keeping above problems in mind. 

The possible uncertainty in cluster age and the MF 
slope, which can be introduced by the data inaccuracy can 
be estimated from the following considerations. The most 
important for our LF-based analysis is accuracy in absolute 
magnitude My . It depends on error in a distance modulus, 
reddening correction, and on a correction for metallicity ef- 
fect. Since no metallicity measurements are published for 
clusters under study we can estimate this correction only 
indirectly. As the cluster galactocentric distances span over 
a range of 7-11 kpc (with the solar value assumed to be 8.5 
kpc), the metallicity radial gradient in t he galactic disk is 
absent according to data on Cepheids of lAndrievskv et alJ 
(2002). Thus we can accept solar heavy element abundance 
and avoid metallicity corrections. The spread of cluster 
metallicity can only be b e due to galactic disk inhomo- 
geneity, which according to IVereschaein fc Piskunovl (Il992f) 
is about A[Fe/H] ~ 0.1. The corresponding absolute mag- 
nitude uncertainty due to metallicity effect is of order of 
AMv ~ 0.1 mag. The standard accu racy of cluster apparent 
dista nce modulus is about 0.2 mag iSubramaniam fc Sagarl 
1999). The reddening uncertainty according to original pa- 
pers is between 0.03 and 0.06 mag and only for NGC 7510 
it exceeds (ue(b-v) = 0.09 mag) the upper limit. These 
cause the My uncertainty to be typically about 0.1 - 0.2 
mag and less than 0.3 mag for the complete sample. Thus 
we believe that aggregated uncertainty in absolute magni- 
tude due to metallicity, interstellar extinction and distance 
modulus errors does not exceed 0.5 mag. This corresponds 
to a random error in logt of the order of 0.15 both for turn 
o ff and tu r n on dating methods. According to simulations 
of lBelikovl ( 19991) one can expect an error in the MF slope 
caused by a My =0.5 mag is of order of 0.2. 

The original CCD photometry was complemented with 
photoelectric magnitudes of brighter stars from the sources 
referred in original papers and/or from the WEBDA 
database. This extends the range of data completeness from 
a faint cluster-specific limit, discussed in Sec. 14.31 to the 
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Figure 2. Maps of stellar fields in the vicinity of studied clusters are constructed with the help of ASCC-2.5, UCAC1, and FONAC 
catalogues. Size of the stellar images corresponds to stellar magnitude (Rfonac f° r th e case of NGC 7510 and R u for the others), 
magnitude scale is shown in the upper right corner. Crosses mark cluster centers. Large circles denote core and corona areas, solid 
quadrangles show used CCD-frames. 



brightest in the observed field stars. Since the data in the 
^-magnitude are the most complete we use in this study 
LFs constructed in V only. 

As one can see from Table Q stellar fields cov- 
ered by the fra mes are comp arable with the cluster 
sizes provided in ILvngal il987f) ca talogue. It should be 
noted however, that ILvngal (Il987l) provides only sizes 
of central parts of open cluster s , wh ich according to 
iKharchenko. Pakulvak fc Piskunovl j200ot) are 2 to 3 times 
lower tha n their full size. Simila r conclusions have also been 
drawn bv lNilakshi et all i2002f) based on the spatial struc- 
ture study of the open star clusters. Thus we should con- 
sider our samples only as inn er population of the clu sters, 
enhanced with massive stars (Kh archenko et aP feoOS). and 
having flatter lu minosity functions and m ass spectra than 
outer population ide Griis et al.ll2002allblci) . 

In Fig.|!2]we show wide neighborhoods of the clusters un- 
der study. They are constructed with t he help of all-sky cat- 
alogues ASCC-2.5 2 llKharchenkcl200ll) for bright est (V < 13 
mag) stars, and UCAC1 3 iZacharias et al"ll2000l) [FONAC 4 

2 available at|ftpj77cdsarc^^trasbg|fr7p^b7cats/I/280A 

3 available at ftp://cdsarc.u-strasbg.fr/pub/cats/I/268 

4 available at ftp://cdsarc.u-strasbg.fr/pub/cats/I/261 



jKislvuk et al]ll999l) in NGC 7510 case] for fainter stars. 
In order to give a reader an idea on a relation of observed 
frames to cluster geometry we show in the maps, positions of 
cluster centers and borders of clu s ter co res and coronae, as 
determined by IKharchenko et alJ <l2003l) from data on spa- 
tial distribution of cluster proper motion members from the 
ASCC-2.5. For NGC 2383 IKharchenko et alJ feOOSt failed to 
detect the core radius and the core is not shown in the clus- 
ter map. Since brightest stars of Hogg 15 are fainter than 
limiting magnitude of the ASCC-2.5, the data on the clus- 
ter cen ter and radius wer e taken from lDias. Lepine fc Alessl 
J2002ft and lMoffaj \l974) respectively. 



One can see from Fig. [5] that regions of two clusters 
(NGC 2383 and Hogg 15) are overlapped with their neigh- 
boring clusters. This should be regarded as a potential 
source of cluster sample contamination and should be kept 
in mind in further discussions. NGC 4755 is the only cluster, 
where observations are carried out in two separated frames 
located on both sides of the cluster center. 
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4 OBSERVED AND THEORETICAL 
LUMINOSITY FUNCTIONS 

Before comparing the empirical and theoretical LFs, it is 
necessary to correct observed LF for different biases men- 
tioned earlier. Among the most serious effects influencing 
observed LFs are contamination of a cluster sample with 
non-members and data incompleteness. Generally the non- 
members are field stars projected onto the cluster area, but 
in some cases there would be stars from a neighboring cluster 
also. The data incompleteness has a strong dependence on 
the crowding/density of stars in the field and the stellar mag- 
nitude such that the data is less complete i n crowded regions 
and towards fainter magnitudes (see e.g. ISaear fc Richtlen 
Il99ll) . All these, as well as details of procedure for the con- 
struction of theoretical LF and fitting of theoretical and em- 
pirical LFs will be described below. 



4.1 Field star de-contamination with galactic star 
count model 

Depending on cluster distance and on the properties of ob- 
served sample, one could apply different methods for the 
selection of cluster members. For nearby clusters, the kine- 
matic method is the most common and reliable. For remote 
clusters the statistical selection with the help of adjacent 
stellar fields are used. Since there is no blank fields observed 
in the neighborhood of the most clusters under study, we 
use the Galactic star count models to estimate the field star 
contamination in their directions. Due to various perturba- 
tions (e.g. due to spiral arms) along line of sight, which can 
be hardly taken into account in details by the model, this 
approach appears to be crude at the first glance, especially 
for a study of the fine structure of the LF. A careful anal- 
ysis indicates contrary to this apprehension however, ft is 
well known that the leading role among a number effects 
influencing the predicted apparent stellar density, is played 
by interstellar extinction, as it is uneven over the galactic 
disk. For example the patchy behavior of the extinction is 
the major reason enabling observation of remote clusters in 
transparency windows. This is a basic point of our approach: 
since the extinction in a cluster direction could be derived 
with much more confidence when compared to an average 
extinction value in the disk, more realistic model predictions 
can be made in a cluster field than in an arbitrary stellar 
field. Below we describe the model components, discuss the 
degree of the model stability against the extinction and stel- 
lar density fluctuations, and compare the model predictions 
with data on selected test fields. 

For each cluster, we compute a theoretical distribu- 
tion of field stars with apparent magnitude ip(V), which 
is specific to the cluster parameters, including galactic co- 
ordinates, interstellar extinction towards the cluster, and 
normalized to the are a of observed cluster region. W e use 
the model described bv|Kharchenko fc Schilbachl<ll996T) . and 
iKharchenko et all (Il997h . The model treats the Galaxy as 
a system, which is symmetric both with respect to rota- 
tion axis, and to the equatorial plane, and consists of three 
populations (thin disk, the disk hereafter, thick disk and a 
spheroid). The population of thin disk consists of several 
sub-populations (MS stars, disk red giants, and super giants 
- stars of I and II luminosity classes) . Each population group 



Table 2. Cluster-specific model parameters 



Cluster 


1° 


b° 


ay 








(mag/kpc) 


Hogg 15 


302.0 


-0.2 


1.00 


NGC 2383 


235.3 


-2.4 


0.38 


NGC 2384 


235.4 


-2.4 


0.41 


NGC 4103 


297.6 


+1.2 


0.59 


NGC 4755 


303.2 


+2.5 


0.92 


NGC 7510 


111.0 


+0.1 


1.13 



is described by its own stellar LF and spatial distribution. It 
is accepted that population LF does not depend on position 
within the Galaxy and is equal to its solar neighborhood 
counterpart. Densi ty distribution s of different populations 
are represented bv iBahcalll II1986T) expressions for disk-like 
and spheroidal subsystems. 

The L F for MS star s bri ghter than My = 13 mag is 
taken from IScalol {l986) and iMurrav et all Jl997l) . while 
for fainter stars (My 1 9 mag) we used data given 
by Ijahreifi fc Wielenl il997T> for nearby stars. The frac- 
tion of evolved stars, which have left MS phase (red gi- 
ants and super g iants) as a function of My was taken 
from IScalol (1986). These values were modified at My ~ 
0.75 mag, (B — V) ~ 1.0 mag acco rding to the data of 
iHolmberg. Flinn fc Lindengrenl lll997l) taking into account 
red giant clump, an equivalent of Population II horizon- 
tal branch for disk-population core helium burning stars. 
The LFs of the thick disk and spheroid and the star num- 
ber density normalization {th in disk:thick dis k: spheroid} are 
{1:0.02:0. 00125} taken fro m iGilmorel lll984l) . Additionally, 
results bv lDa Costal <ll982t) for globular clusters were used to 
model the influence of the horizontal branch on the spheroid 
LF. 

We adopt the following general parameters of the 
model. The galacto-centric distance of the Sun is equal to 
8.5 kpc; the ratio of axes of the spatial density distribution 
for the spheroid is equal to 0.85; the length scale of the disk 
subsystems and the height scale of the thick disk are equal 
to 4 kpc and 1 .3 k pc respective ly. We assume according to 
ISchmidtl lll963ft and lScald lh986Tl that the height scale of the 
disk for MS stars rises from 90 pc to 350 pc at My interval 
[2.3, 5.4] mag, whereas for red giants, it changes from 250 pc 
to 400 pc at My = [-0.75, 2.6] mag. The height scale of 90 
pc was used for extremely young population of super giants. 

Young open clusters are located in the galactic disk with 
irregular interstellar absorption and their extinction param- 
eters differ from average Galactic values. This is a major rea- 
son for using in model calculations a clust er-specific absorp - 
tion values ay, which were derived from IPaxenagd <ll940T) 
formula 



Ay = 



ay 



j sin b\ 



exp 



-d- 



ib\ 



where d is cluster distance from the Sun, b is cluster galactic 
latitude, Ay = 3.1 E(B - V) with E(B - V) taken from Ta- 
ble and hz is height scale of the extinction layer assumed 
to be equal to 100 pc 5 . Values of ay together with cluster 
galactic coordinates are shown in Table H 



u bee paper of I Kilpio fc Malkovl (1997) for recent discussion of 
the extinction formula parameters. 
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Figure 3. A comparison of stellar counts in test fields, neigh- 
boring to NGC 2383/4 (a), and NGC 7654 (b) clusters, with a 
prediction of our model. The histograms show empirical data, the 
bars are the Poisson deviations, indicating data uncertainty. The 
curves are the model counts. 

In order to illustrate the decisive role of the reddening 
effect in the model counts construction, we have confronted 
our results to calculations with a standard value of galactic 
specific extinction a v = 1.6 mag/kpc, adopted for the galac- 
tic disk. As shown by the comparison, for all clusters in the 
working magnitude range of V — 10 — 20 mag, the overes- 
timation of specific extinction leads to considerable under- 
estimation of the model density. The corresponding model 
curves differ in working magnitude range for different clus- 
ters by A log N = 0.3 - 0.8 (i.e. by 70 to 180 %). The model 
stability with respect of observed errors in the value of spe- 
cific extinction was evaluated by variation of values of av 
from Table[3]by 0.2 mag/kpc. According to Sec.[3]this value 
can be regarded as upper limit of error in specific extinction 
values in our sample. We have found that in this case the 
model counts variations are of order of |AlogiV| < 0.15 for 
all clusters under study. 

In order to estimate the effect of spatial density fluc- 
tuations, we consider galactic spiral arms as the most 
pronounced galactic disk sub-structures. We have carried 
out model calculations for cluster areas with and without 
the effects due to spiral arms. We considered three local 
arms (Sagittarius-Carina, Orion and Perseus ones) repre- 
sented by logarithmic spiral s with parameters taken from 
iMarochnik fc Suchkovl l| 19841) . It was assumed, that the den- 
sity input from the spiral structure is equal to yD(R), where 
D(R) is regular density profile in the disk, and 7 is a scaling 
factor, taken to be equal to 0.1, 0.2, and 0.3. The model 



calculations show that for all clusters in the working mag- 
nitude range, the influence of the spiral arms on the stellar 
density is negligible. The corresponding increase of AlogAf 
is about 0.05 for 7 = 0.3, and less than 0.02 for 7 = 0.1. 
These values indicate that the effect of density fluctuation 
is much weaker than the extinction. 

At last we compare our model calculations with ob- 
served stellar counts in selected stellar fields. Unfortunately 
among the clusters of our sample only observations of 
NGC 2383/4 are supplied with photometry in the blank field 
containing about 300 stars. The field is about 5 'away from 
the clusters and thus is partly overlapped with NGC 2383/4 
areas (cf. Fig. 01. For model calculations we used data for 
this cluster from Table |5] Corresponding distributions are 
compared in Fig. I^a). One can observe reasonable quali- 
tative agreement of empirical and model distributions. Ac- 
cording to x 2 -test the both distributions agree at 80-percent 
significance level. Small and generally insignificant access of 
brighter stars (V < 16 mag) can be attributed to the ad- 
mixture of NGC 2384 stars. 

An excellent opportunity for verifying of our model pro- 
vides wide-field CCD-photometry of more than 17 000 star s 
in the field of NGC 7654, published bv lPandev et alJ <200ll) . 
The photometry covers an area of 40x40 arcmin 2 , centered 
at the cluster with a size of about 24'. This leaves a wide 
circular blank field area, enabling not only to confront the 
model counts and observations, but also to estimate the 
stochastic variations of circumcluster fields themselves. It 
is worth to note, that NGC 7654 is located only about 1 
degree North of NGC 7510, a cluster of present sample. 

In Fig. CUb) we show observed and m odel distribu- 
tions f or NGC 7654 blank field area. Although lPandev et alJ 
(2001) do not provide explicit value of the cluster radius, 
their density profiles show that there is no sufficient clus- 
ter population at ro > 13.'3. For construction of the ob- 
served distribution we use all stars located at r > ro, but 
the stars from the SE segment of the circumcluster area are 
not considered. This is where the neighboring cluster Cz- 
ernik 43 having radius of about of 10 'is located, according to 
iKharchenko fc Piskunovl (1200311 . The cluster-specific model 
parameters (d = 1380 pc, and ay = 1.80 mag/kpc ) were 
taken from data for NGC 7654 o flPandev et al.l 1I2OOII) . Since 
there are indications of outward increase of t he cluster red- 
dening , we used the upper limit provided bv lPandev et alJ 
(2001) values. As one can see from Fig. |3fb) the agreement 
between the model and observations is good. According to 
X 2 -test the both distributions agree at 90-percent signifi- 
cance level in the range V =12-19.5 mag. In order to es- 
timate the degree of stochastic variations of circumcluster 
population we arbitrarily divided the whole blank field into 
two equal parts and compared them to the united field sep- 
arately. The x 2_ test has shown that their statistics is some- 
what different: they are the same as the total distribution 
on 85-percent level. This indicates that the model counts 
approach to the member selection is at least not worse than 
the blank-field method based on an arbitrary sampled test 
area. 

Theoretical differential distribution of stars with appar- 
ent magnitude ip(V), derived with help of this model were 
used for construction of reduced LFs as it is described in 
Sec. Em 
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Table 3. Contributions of field star and neighboring clusters 
population are compared with the projected stellar densities of 
NGC 2383 and Hogg 15 clusters 



Population 



Projected density (l/d°) Designation 





NGC 2383 


Hogg 15 




Observed (total) 


2724 


6077 


^cluster 


Field stars 


932 


1870 




Ncighbor+field 


1273 


2401 


^neighbour 


Neighbor 


341 


531 


^neighbour 


Cluster 


1451 


3676 


^cluster 



4.2 Contamination by an overlapping cluster 

The contamination by an overlapping cluster is studied using 
the projected density derived from stellar counts in recent 
all sky catalogues. The star counts in the USNO-A2 6 cata- 
logue is strongly affected by crowding effect in the areas of 
Hogg 15 and NGC 2383/2384 clusters (which is not surpris- 
ing for Schmidt camera based surveys in dense fields). To 
avoid crowding effect we used astrograph based catalogue 
UCAC1. This catalogue is a part of all-sky survey UCAC, 
covers Southern Hemisphere south to feooo ~ —20°, and 
contains the fields under study. 

NGC 2383/2384 lie at the very edge of the survey 
UCAC1 area, where only short exposures provide star counts 
homogeneity and completeness at R u < 13.7 mag. We con- 
structed density profiles in a field centered at the NGC 2384. 
In order to estimate the proper stellar density in a cluster, 
we counted stars in four quadrants separately. The average 
observed star densities in the field E/ and in NGC 2384 
corona areas E2384 were determined from the counts in the 
first and forth quadrants of the corresponding areas located 
at 02000 ^ 7 h 42 and within 0?10 -f- 0?15 and 0?35 -5- 0?60 
from the cluster center which are free from contamination 
of NGC 2383 members. As the values of E/ and E2384 do not 
strongly vary over the studied area, the density of NGC 2383 
cluster members can be calculated as 02383 = E2383 — E2384 
(E2383 is counted within NGC 2383 frame). The star density 
of the corona of NGC 2384 itself, then is 02384 = E2384 — E/. 
The results are shown in Table [3] and one can see that 
contamination of NGC 2383 member sample with stars of 
NGC 2384 cluster computed as <T2384/o"2383 is about 24%. 

We cannot estimate the contamination of NGC 2384 
member sample with stars of NGC 2383 cluster as there are 
no stars with R u < 13.7 mag within NGC 2384 frame piece 
overlapping with NGC 2383 cluster. Moreover, this frame 
piece is small and we can assume insignificant number of 
NGC 2383 members there. 

Similar procedure was applied to Hogg 15 cluster. The 
counts were performed down to R u — 15 mag. The average 
observed density in corona of NGC 4609, E4609 was com- 
puted from the stellar counts at distance 0?15 — 0?20, and 
average observed density of the field stars at 0?35 — 0?60 
from the center of NGC 4609. As one can see from Table |3] 
Hogg 15 is contaminated with NGC 4609 stars by about 
14%. This fraction is smaller than in the NGC 2383 case, 
but is not negligible. 



6 Monet D., Bird A., Canzian B. et al., USNO-A 2.0:A Catalog of 
Astrometric Standards (CD-ROM distribution). US Naval Obs., 
Washington, 1998. 



4.3 Incompleteness effects 

It is well known, that completeness of a photometric sur- 
vey is achieved only at about 2 mag brighter than Vu m (so 
called completeness limit V cmp ). Below V cmp the observed 
LF does not reflect behavior of a genuine LF, and this mag- 
nitude range should be excluded from the consideration. It 
is also important to be confident, that for clusters, having 
observations in two CCD frames both fields have the same 
completeness limit. Unfortunately, no V C m P determinations 
were made in original papers. This forced us to determine 
the Vcmp within this study. Since no test fields were ob- 
served in original papers we are only able to estimate roughly 
the completeness limit as a position of apparent maximum 
in the distribution of stars with V magnitude within ob- 
served frames. Since the incompleteness increases gradually, 
the above estimate should be considered as upper limit of 

The estimated values of V C m P are shown in Tabled For 
NGC 4755 we find good agreement between V cm p in both 
observed frames. Note, however, that the frames do not in- 
clude the central area of the cluster, containing the brightest 
cluster members. This causes another kind of incompleteness 
i.e. incompleteness for bright stars. One should keep this in 
mind when the final result for this cluster is discussed. In 
contrast the frames of other clusters are located in central 
regions and no bright star incompleteness is expected there. 

The other kind of incompleteness, which should be con- 
sidered is incompleteness due to overlapping of stellar images 
in the CCD frame, which we call here as "c rowding effect 
incompleteness" (see ISaear fc Richtlerlll99ll for details). It 
arises due to finite size of stellar images which depend on 
stellar brightness. An image of a star of some brightness oc- 
cupies certain area within the CCD frame, hides some neigh- 
boring stars of lower brightness, and excludes them from the 
statistics, affecting the brightness function ip(V), a distribu- 
tion of stars with apparent magnitude. The relation between 
corrected for the crowding effect brightness function ij) ce (V) 
and observed one ip a (V) is expressed as 



ip ce (V) = MV) 



1 + 



s(v')Mv')dV 



Here V m in is the magnitude of the brightest star in the 
frame, s(V) is an area occupied in the frame by an image of 
a star of magnitude V. Both s(V) and ^> C e(V) are expressed 
in units of the frame area. We have evaluated the crowding 
effect for each cluster of our sample. The functions s(V) were 
found empirically for every frame from the image statistics 
as s(V) = Tvrp, where r p is average minimum distance in 
the plane of a frame from stars of magnitude V to remaining 
stars. As we found for the considered range of magnitudes 
this effect only negligibly (less than two percents) changes 
distribution ip a {V). Thus we assume ip C e(V) = ip a (V) in our 
further discussions. 



4.4 Theoretical luminosity functions 

For construction of a theoretical LF <f>(Mv), we have as- 
sumed a model of continuous stars formation. It is computed 
as 
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Table 4. Open cluster parameters derived from the fit of the LFs. 



ct>(M v )= [ MMv)\(t)dt, 

J t 

where to and t\ are minimum and maximum ages of the 
cluster stars, function 4>t(My), computed using equation 
is the LF of stars with age t (to ^ t ^ ti), X(t) is the star 
formation rate (SFR) at age t. The value of t\, representing 
duration since formation of first stars, could be regarded as 
the cluster age. 

The mass - absolute magnitude relation m(Mv,t) and 
its derivative were calculated along the isochrone of age t 
using a cubic spline interpolation. For the IMF we considered 
a power- law representation f(m) = km~ x , with k to be 
normalizing factor, and x the IMF slope to be determined 
within this study. We assumed a constant SFR in our model, 
X(t) = const. The resulting parameters for each cluster were 
drawn from the best fit of the theoretical and observed LFs. 

In order to construct theoretical isochrones and LFs 
which include both Post-MS and Pre-MS stages for ages 
typical to clusters of our sam ple, we combined Population I 
Pre-MS evolutionary tracks of iD'Antona fc Maziteul (Il994f) 
for masses 0.1 to 0.8 m®, and of lPalla fc Sta hlcr ( 1993|) for 
masse s 0.8 to 6 mp> and Maeder' group Post-MS calcula- 
tions dSchaller et al.lll992l) for m = 0.8 - 120 m©. The grids 
were properly tuned to provide a continuous transition from 
Pre- to Post-MS ages and smooth and uniform mass - lu- 
minosity and mass - radius relations along the ZAMS. The 
isochrones were computed from the models corresponding 
to the Population I chemical abundance (Y, Z) =(0.30,0.02) 
using linear interpolation. 

In order to convert the theoretical isochrones from 
log T e ff, log L/Lq plane to the observed (B—V)o, My plane, 
we used bol ometric corrections an d (B — V)o — logT e // re- 
lations from ISchmidt-Kalerl Jl982l) tables for the luminosity 
classes I, III and V. 



4.5 Fitting of theoretical and observed LFs 

The following iterative steps are used in the fitting proce- 
dure: 

(i) Apparent LF construction: Apparent LF (observed 
brightness distribution) ip a was constructed from the data 
available for given frame(s) as smoothed density estimation 
with a rectangular one-magnitude wide window, and a step 



A log* 

0.0 ±0.2 

0.2 ±0.4 
0.2 ±0.4 
0.2 ±0.4 
0.6 ±0.6 
0.2 ±0.6 
0.2 ±0.2 
0.5 ±0.7 
the upper MS isochrone fitting 



of 0.25 mag. The histogram form of smoothing kernel was 
selected to apply correction for field star contamination. 

(ii) Construction of a cluster LF: We assume that appar- 
ent LF ip a (V) is a composition of a distribution of cluster 
stars with apparent magnitude tl> c (V) and field stars bright- 
ness function represented by galactic star count model dis- 
tribution ipf(V): 

where free parameter p denotes percentage of field stars, 
contaminating the sample. With help of apparent distance 
modulus V — Mv the function ip c is transformed to absolute 
magnitude scale, and is called hereafter as observed cluster 
LF <t> D {Mv). 

(iii) Luminosity function fitting: Theoretical LF <f>(My) 
smoothed in the same way as observed LF was fitted to 
empirical function <f> c (My) within a cluster-specific range of 
magnitudes. The entropy 

A = \MMy) log(^ c (A4)Mil4))| 

was constructed to find the best-fit parameters {p, x,to,ti}. 

Steps (ii) and (iii) were repeated iteratively by varying 
p until the best agreement (in terms of \ 2 statistics) between 
theoretical and observed LFs is achieved. The parameter p 
was varied between and 1 in a step of 0.1. 

Internal accuracy of a;,logt = logii, and Alogi = 
logti — log to was estimated on the basis of the ke rnel 
smoothing theory dSilvermanl Il986l . lLapko et alj Il996l) for 
selected histogram grid parameters (range and step). An 
accuracy of the field stars percentage p was estimated from 

2 _ E(p»-f) 2 /x 2 -_Ep 1 M_ 
EVx? ' P EVxf 

assuming that o~p is a second order central moment of a 
distribution function equal to that of 1/x 2 - 



5 DISCUSSION 

The results of the LF fitting are listed in Table [I] For each 
cluster we show fit parameters: the range of absolute magni- 
tudes selected for the LFs fit, the best-fit \ 2 value achieved 
in the iterations, and derived parameters: position of H- 
maximum in theoretical LF My , indicating absolute mag- 



Cluster 


Fit parameters 




Derived parameters 




My range 


x 2 


My 


P 


X 


logt 




(mag) 




(mag) 








NGC 2383 


-6.0, 3.0 


23.6 


2.50 


0.1 ± 0.2 


2.2 ± 0.4 


7.1 ±0.2 


NGC 2383 a 


-2.0, 3.0 


37.2 




0.2 ± 0.2 


1.2 ± 0.4 


8.5 


NGC 2384 


-5.0, 3.0 


39.6 


2.00 


0.3 ±0.2 


1.0 ±0.5 


7.0 ±0.2 


NGC 2384 6 


-5.0, 3.0 


41.5 


2.00 


0.2 ± 0.2 


1.0 ±0.5 


7.0 ±0.2 


NGC 4103 


-3.5, 2.0 


24.1 


1.25 


0.3 ± 0.2 


1.5 ± 0.6 


6.8 ±0.1 


Hogg 15 


-5.0, 2.5 


79.2 


0.75 


0.6 ± 0.2 


1.3 ± 0.4 


6.7±0.3 


Hogg 15 c 


-5.0, 2.5 


17.0 


0.75 


0.5 ± 0.2 


1.3 ± 0.4 


6.6 ±0.3 


NGC 4755 


-3.0, 3.0 


10.8 


2.50 


0.2 ± 0.2 


1.4 ±0.3 


7.2 ±0.2 


NGC 7510 


-6.0, 1.5 


13.8 


0.75 


0.3 ± 0.2 


1.2 ± 0.4 


6.8 ±0.3 



a corrected for contamination by NGC 2384 stars, logt = 8.5 is accepted from 
b derived in the frame area free from NGC 2383 stars 
c corrected for contamination by NGC 4609 stars 
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nitude of the turn on point, derived percentage of field star 
contamination p, the IMF slope x, adopted age of a cluster 
logt, and age spread parameter A logt. Note, that parame- 
ter errors shown in Table 2] are in good agreement with the 
uncertainty estimations given in Sec. [3] on the basis of data 
accuracy analysis. The results are discussed below separately 
for clusters under study. 

In Fig. 0] we show the CMDs of clusters under dis- 
cussion. The reddening and distance parameters are taken 
from Table Q and fine tuned by the variation of the tabu- 
lar values within their accuracy to reach the best agreement 
of the cluster CMD, empirical ZAMS, and corresponding 
isochrones. Cluster proper motion members and non mem- 
bers are marked, if cluster stars are bright enough to be 
included in the ASCC-2.5. To give an idea how the adopted 
distances and reddening values agree with the photometry 
we show in th e diag ram a position of empirical ZAMS of 
ISchmidt-Kalerl il982T) . In order to illustrate how ages de- 
rived from the LF analyses confirm with cluster CMDs, we 
show in Fig. 0] the isochrones corresponding to the ages, to 
and t\. 

In Fig. we display the LFs. Observed cluster LF <j> c , 
constructed as described in Sec. l4.5l is shown with filled his- 
togram, and corresponding theoretical LF <f) fitted to (f> c is 
shown with a curve (the fitted portion of (j> is shown with 
heavy curve, while the rest is marked with thin curve). For 
comparison purposes we show also apparent LFs ip a (thin 
histogram and hatched area), and contribution from field 
stars tpf (dotted curve), both displaced by the value of 
V — My- Below we discuss each cluster as per their ap- 
pearance in figures 2] and 



5.1 Hogg 15 

In the cluster CMD, one can indicate the MS turn on point 
at M v ~ 0.5 mag. Note, that the H-feature should be also 
found in the LF in the vicinity of M v , where apparent LF 
of Hogg 15 stars also shows a local maximum. The data 
incompleteness dominates after My = 3 mag, where ip a 
gradually decreases in agreement with Table ^ indicating 
that M v mv — 3.1 mag. The field star contamination for 
Hogg 15 is found to be highest amongst clusters considered 
here. This effect, as it is seen from the filled histogram in 
Fig keeps the LF practically unchanged at the brighter 
portion (Mv < 0.5 mag), and hides substantially the H- 
feature region. 

In order to take into account contamination from an 
overlapping cluster NGC 4609 theoretical LF <f)(Mv) was 
composed of two populations belonging to these clusters. 
Since the clusters reside at different distances from the Sun, 
the stars of NGC 4609 should be shifted with respect to 
that of Hogg 15 by the difference of cluster apparent dis- 
tance moduli A(V - M v ). Since V - M v of the NGC 4609 
according to LGM2.2 is equal to 11.45 mag the value of 
A(V — My) is taken to be equal to —5.05 mag. Due to small 
fraction of NGC 4609 stars projected on the Hogg 15 area 
(about 15% of Hogg 15 population according to Sec. 14. 21 
the specific shape of the corresponding mass spectrum is 
not important and we assume that it follows Salpeter l aw. 
The age of NGC 4609 was taken from iMermilliodl (ll98lT) as 
logt = 7.56. Composite LF is shown in Fig. with solid 



curve, while the input from the NGC 4609 population is 
shown with the dot-dashed curve. 

As one can see from Table 2] taking into account pop- 
ulation of NGC 4609 does not strongly change cluster pa- 
rameters (except the age spread), but considerably improves 
agreement between theoretical and observed parameters no- 
tably reducing \ 2 parameter. 

Fig. 2] shows that ages <o,ii derived from the LF fit- 
ting are in good agreement with cluster CMD. The observed 
cluster LF 4>c also demonstrates a high degree of conformity 
with composite theoretical LF. A deep and narrow gap of 
(f> c at about Mv ~ 2 mag shows that the derived input of 
field stars p is somewhat overestimated. As shown in Ta- 
ble [I] Hogg 15 is one of the two clusters in our sample, 
where considerable (although quite uncertain) age spread of 
A logt = 0.6 is detected. This value, however, is reduced to 
the insignificant value of A logt = 0.2, which is common to 
the whole sample, if one takes into account the contamina- 
tion from NGC 4609. We therefore conclude that the age 
spread in Hogg 15 is introduced mainly due to admixture 
of field and overlapping neighbor' stars, and consider the 
cluster parameters as derived for the Hogg 15 c case. 

Hogg 15 is the youngest cluster of our sample, with an 
age of 4 Myr. How does this value compare with the earlier 
published values? For the first time, the age of Hogg 15 
was es t imate d by iMoffad lll974h as 8 Myr. According to 
iLvngal lll987l) the age is 10 Myr. Recent values come gen- 
erally from isochr one fitting o f the u pper MS: LGM2.2 de- 
fine it as 6 M vr. ISa ear et al] j200l|) provide age estimate 
of 6 ± 2 Mvr. lPiatti fc Clarial ihoOll) . who found initially 
it as 300 My r, later reduced t his unusually high value to 
20 ± 10 Myr JPiatti et al.ll2002h . Technique of matching in- 
tegrated theoret i cal an d observed spectra of the cluster lead 
lAhumada et alJ j200Cf) to the age values 5 ± 2 Myr or 30 
Myr. Since MS turn off ages are biased toward higher val- 
ues, our age estimate, based on MS turn on point technique, 
is in fact the lowest one. Note, however, that the majority 
of the above mentioned values agree with the present one 
within their acc uracy. The IMF sl ope was determined for 
Hogg 15 stars bv lSagar et alJ B2001T) . who find it to be equal 
1.35 ± 0.21, which again coincides with the present result. 

5.2 NGC 7510 

It can be seen from the CMD that the regular cluster se- 
quences are embedded in a "cloud" of presumably field stars. 
The astrometric members also show definite spread around 
the MS turn off point. It may be that some of them are co- 
moving field stars. The MS turn on point can be detected 
somewhere between M v ~ 0.5 and 1.5 mag. The apparent 
LF ip a shows a broad plateau between My » 0.5 and 3 
mag, with weak maximum at Mv = 2.25 mag. As indicated 
by galactic model star counts, it should be formed mostly 
by field stars. The sample demonstrates certain incomplete- 
ness below the maximum, where ip a instead of steady in- 
crease slowly falls down. According to Table^the complete- 
ness limit of original data is M v mp = 0.9 mag. We select 
Mv = 1.5 mag as faint limit of the LF fitting range based 
on the value of \ 2 parameter. The agreement between ob- 
served cluster LF ip c and theoretical LF is one of the best 
amongst the studied clusters. 

As one can see from Fig. [I] the ages to = 2 Myr and ti = 
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Figure 4. Color-magnitude diagrams of studied clusters. Open circles are stars from original observations corrected for adopted V — My 
and E(B — V) values. S quares and cros s es are proper motion cluster members and field stars respectively from ASCC-2.5. Heavy curve 
is an empirical ZAMS ofjSchmidt-Kalcr 1 1982), and lower and upper thin curves are isochrones corresponding to derived age parameters 
tl and to respectively. For NGC 2383 we show for comparison both isochrones of logt = 8.3 and 8.5, fitted to its upper MS, and 
isochrones of logt = 6.8 and 7.0, corresponding to age parameters of NGC 2384, derived from the analysis of its LF. For Hogg 15, we 
show corresponding isochrones as listed in Table|4] The long dash curve is the MS of NGC 4609 as it would appear for Hogg 15 reddening 
and distance. 
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6 Myr, derived from the LF fitting are in good agreement 
with details of the cluster CMD. The lower pre-MS branch 
of the cluster coincides with ti isochrone. The to isochrone 
fits the apparent sequence of stars deviating from the MS at 
My ~ 1 mag. Stars forming this group do not show evident 
concentration to the cluster center. If, however, this group 
is not a random concentration of field stars, one can regard 
this as an evidence of the second star formation event in the 
cluster. 

Excellent agreement between theoretical and observed 
LFs concerns only MS. We are unable to follow the H-feature 
over its full width, and are not certain if the descending por- 
tion of the ipc is related to this detail or it is a consequence 
of the data incompleteness of the present survey. Note, how- 
ever, that cluster age derived from the assumption that H- 
maximum is located at My ~ 1 mag (just coinciding with 
the data completeness limit) is in agreement with the CMD 
arguments, which is independent of data incompleteness. In 
contrast, the mass spectrum slope, found from brighter por- 
tion of the LF should be regarded as a confident one . 

MS turn off age of NGC 7510 can be found in iLvngH 

<1987t) catalogue ( 1 My r), and in LGM2.2 list (38±2 Myr), 
ISaear fc Griffiths! <199lb find the age < 10 Myr. Again 
our age of 6 Myr is the lower estimate of these values. 
Mass function of stars of NGC 7510 was constructed by 
ISaear fc Griffiths! ( 1998) . It was found, that in the mass 
range 1-14 rriQ the mass function slope a; is 1.1 ± 0.2, which 
again fits present results well. 

5.3 NGC 4755 

NGC 4755 is the only cluster of our sample, which is ob- 
served in two frames, located on both sides of the dense 
cluster center (Fig. We find that both the frames are 
similar with respect to photometric quality and the com- 
pleteness issue. So we consider them together. 

The cluster CMD shows well defined stellar sequences. 
The MS turn on point can be seen at My « 3 mag, and cor- 
respondingly H-feature is expected to be located at My — 3 
mag and fainter. The uncorrected LF i/j a shows a step at 
My — 2 — 3 mag, related to cluster H-feature. The data 
incompleteness dominates after My = 5.5 (5.2 according to 
Table0 mag, where i/> a gradually falls down. Correction for 
field stars for this cluster can be made with a certain confi- 
dence. By comparing with model brightness function ipf, it 
can be seen that at brighter magnitudes the incompleteness 
can be regarded as negligible. Thus selecting My = 3 mag 
as a faint limit for the LF fitting range, we are safe from 
the data incompleteness bias. Since the brightest stars, re- 
siding in the cluster center and shown in Fig. Q>] are not 
present in original data (see Fig. |3J, we select My = —3 
mag as the bright limit of fitting range. As it is seen from 
the filled histogram of Fig|K| we correctly identify the H- 
maximum position, and in spite of the fact that the fitting 
was made only for MS portion of the LFs, the theoretical 
LF reproduces well the pre-MS H-feature after LF reaches 
its minimum and turns over again. However, we are not able 
to reach the convective portion of the LF due to relatively 
brighter completeness limit of the present data. 

It is seen that theoretical LF is in good agreement with 
the the observed one <f> c even outside the fitting range. The 
derived age spread of cluster stars, according to Table [I] is 



insignificant, and we conclude that no evidence of continuous 
star formation can be found from the LF analysis for this 
cluster. As one can see from Fig. 0] the age of 16 Myr as 
derived from the LF fitting technique is in good agreement 
with the cluster CMD. The IMF slope coincides within its 
accuracy with the Salpeter value. 

MS turn off age of NGC 4755 can be found in iLvngal 

Jl987l) catalogue (7 or 24 Myr ), and in LGM2.2 list (16 ± 3 
Mvrl. lSagar fc Cannon] dl995l) find that it is about 10 Myr, 
and most of pre-MS stars have ages between 3 and 10 Myr. 
Recent age estimation 10 ±5 Myr, derived f rom the fitting of 
upper MS, is given bv lSanner et alj i200 lT) . For this cluster 
our age of 16 Myr is rather average than the lower estimate 
of published values . Mass function of t he NGC 4755 stars 
was constructed bv lSanner et alJ |2001). It was found, that 
in the mass range 1-13 m© the slope x is 1.68 ± 0.14, which 
is significantly steeper than the value estimated here. 

5.4 NGC 4103 

The observed frame (see Fig. |5J covers only a minor part of 
the cluster center. One can therefore expect the results to 
be more uncertain in comparison with the majority of other 
clusters under study. 

In the cluster CMD, the MS turn on point can be seen 
at My ~ 1.75 mag, and the H-feature dip is expected at 
My = 2 — 4 mag, seen as the CMD area of lower density. The 
apparent LF ip a shows week local maximum at My = 1.5 
mag. The data incompleteness (see Table dominates at 
My > 4.8 mag, where ip a gradually falls down. However, 
%p a flattens after My = 4 mag, which indicates existence of 
certain data incompleteness at brighter magnitudes. To be 
completely safe against this effect and taking into account, 
that the shape of ip a at My = 2.5 — 3.5 mag is very close 
to that of model field star brightness function, indicating 
high fraction of field stars in this magnitude range, we se- 
lect My = 2 mag, as a faint limit of the LF fitting range. 
The minimum of % 2 is achieved at p = 0.3, with the corre- 
sponding parameters shown in Table 2] 

In spite of the cluster LF <j) c turning to increase af- 
ter My ~ 3 mag, we could not regard the H-feature with 
same confidence as in Hogg 15 or NGC 4755, due to poor 
statistics. It should be noted, however, that as in the case 
of NGC 4755, theoretical LF, being fitted to the MS obser- 
vations, reproduces the empirical Pre-MS stage reasonably 
well. Cluster age, derived from the LF and MS turn on point 
analysis shows, that this cluster seems to be younger, than 
it was regarded before, generally from MS turn off dating 
technique (see below). It should be noted, however, that 
NGC 4103 is the case, where the MS turn off technique due 
to practically unevolved upper main sequence of the cluster 
is able to provide an upper age estimate only. In contrast the 
mass spectrum slope, found from brighter portion of the LF 
should be regarded as accurate enough. Note some excess of 
brightest stars {My ~ —3 mag), which can be however, a 
consequence of random placement of rare bright field stars 
over the cluster region. 

All the authors provide NGC 4103 MS turn off age 
and are u nanimous in i ts value. The age is 22-41 Myr ac- 
cording tolLvngal lll987l) . 24 ± 9 Myr according to L GM2.2 
ISagar fc Cannon! (Il997ft find it to be about 30 Mvr. lForbesl 
l^96TT "estimated about 25 Myr and notes, that MS turn 
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Figure 5. Luminosity functions of the clusters under study. Histograms show observed LFs: an apparent LF ij) a (My) (hatched), and 
cluster stars LF ipc(My), corrected for the effect of field stars (filled), bars are Poisson deviations, indicating data uncertainties. Curves 
are theoretical LFs: cluster star LF cj>(My) fitted to (f> c (My) (solid), and field star brightness function ip(V) displaced along the abscissa 
by V — My (dashed). Heavy portion of the theoretical LF marks the range of fitting. For Hogg 15 and NGC 2383 solid curves show 
composite LFs, which take into account the influence of a neighbor, and dash-dotted curves show model LF, which mimics neighboring 
cluster star population (sec Sec. I5.1l and l5.5l for details). 



14 A.E.Piskunov et al. 



on age shou l d not be much different, and according to 
ISanner et alJ I^OOl) the age is 20 ± 5 Myr. Again we give 
the lowest value of the cluster age of about 6 Myr and be- 
lieve, that both cluster CMD and LF support this value. 
Mass function of s tars of NGC 4103 was constructed by 
ISanner et all (1200 ll) . It was found that in the mass range 
0.7-12 m,Q the slope x is 1.46 ±0.22. Almost the same value 
is estimated here. 

5.5 NGC 2383 and NGC 2384 

In NGC 2383/2384 pair, the dominating role in the neighbor 
contamination is played by the extended and loose cluster 
NGC 2384 (see Sects. ETTl |4"2"1 . In contrast, the central part 
of NGC 2384 can be regarded as relatively free from the 
NGC 2383 influence. Therefore we start this section with 
discussion of NGC 2384 results. 

The CMD of NGC 2384 in its lower part is very 
fuzzy, while the MS part of the diagram is well defined. 
ISubramaniam fc SaeaJ lll999l) even assume that stars fainter 
than V = 16 (My = 2.8) mag are mainly field stars, as 
can be seen from the Fig. |1] Nevertheless, a certain frac- 
tion of cluster members is observed even at faintest mag- 
nitudes. These are pre-MS stars, forming the upper bound 
of the observed sequence at My > 2 mag. In contrast to 
other clusters the CMD of NGC 2384 does not display an 
extended halo of field stars, where cluster MS- and Pre-MS 
sequences are embedded. Field stars (presumably giants and 
late type dwarfs) form a compact clump in the bottom of 
the diagram, with (B — V)o ~ 0.6 mag and Mv > 4 mag, 
leaving Hertzsprung gap to be unoccupied. The MS turn on 
point in the CMD is less certain than for other clusters. If 
one correlates cluster CMD with the morphology of the ob- 
served LF one should detect it around M v = 2 — 2.5 mag, 
where cluster MS-stars certainly cease to appear, and simul- 
taneously a local detail, resembling the H-feature appears in 
the observed LF. One can observe, however, that a few of 
the Pre-MS stars do appear in the CMD at brighter magni- 
tudes (My ~ 0.5 — 1.5 mag). Their extension forms an upper 
bound of a broad Pre-MS branch, extending down to the 
limit of observations. The regularity of this feature implies, 
that it is not a random pattern of field stars. They would be 
rather unresolved Pre-MS binaries, or the youngest gener- 
ation of NGC 2384 stars. Note, however, that this brighter 
position of the MS turn on point does not produce anything 
resembling the H-maximum, expected to be located near 
this point (see Sec.|5Jl. 

The LF of NGC 2384 demonstrates irregular behavior, 
and its agreement with theoretical LF is worse in compar- 
ison to the other clusters under study. This is confirmed 
statistically as the x 2 value is highest for this cluster (ex- 
cept Hogg 15 case with no correction for contamination by 
NGC 4609). One can mention an excess of observed num- 
ber of massive stars (Mv ~ — 4 mag), and a deficiency of 
medium-mass stars (Mv = — 1 -5- 1 mag) with respect to 
theoretical LF 0. Also LF-detail, which we identify with the 
H-feature is poorly developed; the observed H-dip is nar- 
rower, than the theoretical one. This behavior is in confor- 
mity with mass segregation scenario in the case if the ob- 
served frame covers the central part of the cluster only. The 
de rived slope of the mass function is similar to values found 
bv lde Griis et al.H2002bl) for central parts of rich LMC clus- 



ters NGC 1805 and NGC 1818. The LF irregularity could be 
attributed then to incomplete member statistics. The higher 
extent of the cluster comparing to the observed frame nat- 
urally explains other peculiarities of NGC 2383/4 pair like 
presence of early type stars in the NGC 2383 field and ex- 
istence of underdeveloped H-feature in its LF (see below), 
or weak excess of bright stars in the neighboring blank field 
(see. Sec. 14. H . 

The stability of NGC 2384 parameters against contam- 
inating effect of NGC 2383 was checked using data from the 
part of NGC 2384 frame, located outside NGC 2383 area 
(see Fig. . As one can see from Table 31 there is no dif- 
ference between derived parameters of NGC 2384 in both 
cases. Thus we can regard the influence of NGC 2383 mem- 
ber contamination on the derived parameters of NGC 2384 
to be negligible. 

Color magnitude diagram of NGC 2383 is more pop- 
ulated than that of the neighboring NGC 2384. There is a 
well defined evolved portion of the MS compatible with clus- 
ter age \ogt — 8.5, whic h is in agreement with a n esti mate 
of logt = 8.45 — 8.6 bv ISubramaniam fc Saearl lll999|) and 
differs from the ages logt = 7.4 and 7.17 derived bv lLvngH 
il987fl and LGM2.2 respectively. This contradiction arises 
due to the presence of bright ASCC-2.5 kinematic members 
with ages of a few tens of Myr. The MS turn on point, corre- 
sponding to this younger age can not be easily identified in 
the CMD. Nevertheless the LF indicates presence of the H- 
detail even in the apparent LF ip a . The MS-band at fainter 
magnitudes (Mv > 2 mag) is widened and unlike the other 
clusters studied here, it is not due to field giants extending 
the MS from the side of fainter magnitudes, but it is due 
to bright stars forming the high luminosity envelope of the 
MS. 

In the corrected LF (f> c , the H-feature becomes more 
pronounced, leaving the H-dip however to be somewhat un- 
derdeveloped (more shallow and narrow than in the cor- 
responding theoretical LF). The cluster age derived from 
the LF fitting corresponds to lower age estimation coming 
from the brightest MS stars. The value of \ 2 indicates that 
the LF fit is more reliable, than in the case of NGC 2384. 
Note the unusually steep value of the IMF slope x = 2.2, 
estimated for this cluster. In spite of formally satisfactory 
results of the LF fitting there are several points, indicating 
that NGC 2383 carries signs of a dual nature, which were not 
taken into consideration. For example, evolved upper MS in- 
dicates that NGC 2383 is about 300 Myr old (see Fig. EJ, 
while LF indicates its age close to that of NGC 2384. The 
relatively young age of NGC 2383 is also supported by the 
presence of bright MS kinematic members. Possibly this is 
the main rea son for assigning younger age to the cluster by 
iLvngal il987h and LGM2.2. This confusion can be removed 
with the help of a model of an overlapping neighbor, with 
the only difference from the case of Hogg 15, that in this 
case degree of contamination with NGC 2384 stars is higher 
(see Sec. 14.21 . This is also supported by independent data 
of the ASCC-2.5, which indicates, that these clusters are lo- 
cated close n ot only in the sky, but a lso have similar proper 
motions (see lKharchenko et alJl2003T) . So a probable proper 
motion member of one cluster can also well belong to the 
other, and the brightest kinematic members of NGC 2383 
could be actually projected members of NGC 2384. This 
can also be supported by the fact that numerous MS-stars of 
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NGC 2383 forming upper envelope at My > 2 mag, could be 
actually Pre-MS stars of NGC 2384. Apart from coincidence 
with positions of corresponding isochrones in the CMD this 
idea is also supported by spatial distribution of these stars. 
They tend to be located in NGC 2383 frame within circle of 
its neighbor. The number of these stars is compatible with 
the degree of contamination derived from stellar counts and 
they might be responsible for steeper IMF slope found for 
NGC 2383, when neglecting influence of NGC 2384. 

In order to take into account contamination from the 
overlapping cluster theoretical LF of NGC 2383, 0(My) was 
composed of two populations belonging to both clusters. Un- 
like to the case of Hogg 15, NGC 2383 and NGC 2384 re- 
side at the same distance within the errors (A(V — My) = 

0. 1 — 0.3 mag). So no correction is required for the LF magni- 
tude scale while constructing the composed theoretical LF. 
We assume, that the mass spectrum of NGC 2384 stars does 
not vary with location and the case of NGC 2384^ in TableH 
can be considered as its representative. The corresponding 
values of to and t\ were also taken from the Table 2] Since 
no developed H-feature exists in the LF at t ~300 Myr we 
were forced to use the LF of NGC 2383 in a customary, way 

1. e. for determination of x and p parameters. Cluster age was 
found from the MS turn off point and fine tuned to satisfy 
the x 2 -test. The lowest x 2 -value was achieved at logt = 8.3, 
which is somewhat less than isochronic age of logi = 8.5 
(Fig.0Jl, determined from the Post-MS isochrone fitting. Due 
to above mentioned reasons we regard this disagreement as 
insignificant and show in the Table [I] the turn-off age. The 
best-fit composite LF for the case NGC 2383 a is shown in 
Fig. with solid curve, while the contaminating model LF 
of the NGC 2384 is shown with the dash-dotted curve. 

As one can see from the Table 3] the proposed model 
provides less accurate fit of theoretical LF to observed LF 
than the previous one. On the other hand, it harmonises the 
full observation scope (star counts, cluster CMDs and LFs). 
As a result, with the age determined from the MS turn off 
point calibration, the IMF slope has been reduced to the 
value close to the Salpeter's one with the value of p almost 
unchanged. 

In the light of above discussions, it is not surprising that 
the cluster parameters derived in the present study are in 
good agreement with those de rived in the original paper by 
ISubramaniam fc Sagarl jl999T) . For example, our values of 
x prac tically coincide with those of ISubramaniam fc Saearl 
(1999), who found x = 1.3 ± 0.15 for NGC 2383 and 
x = 1.0 ± 0.15 for NGC 2384 Similarly, in agreement with 
ISubramaniam fc Saga] (Il999l) . we regard that NGC 2383 
ages listed bv lLvngal dl987l) as 25-41 Myr or by LGM2.2 as 
16 Myr are underest imated, whil e the NGC 2384 ages (1- 
10 Myr according to lLvngalll987i and 8 ± 1 Myr according 
to LGM2.2) are in fair agreement with our result. However, 
again the age of NG C 2384 derived from LF-fittin g is lower 
than that derived bv lSubramaniam fc Sagarl il999T) from the 
upp er part of the MS. 

ISubramaniam fc Sagarl 1^999) conclude that in spite 
of their spatial proximity clusters NGC 2383/2384 do not 
constitute a physical pair in the sense of common origin. 
Present study not only supports this point but makes it more 
stronger. In spite of the fact that both clusters do overlap 
in the plane of sky, reside approximately at the same dis- 
tance from the Sun, and have similar proper motions, they 



are very different in all other respects e.g. cluster geometry, 
morphology, environment, stellar contents, and age and can 
not be regarded as twins. 



6 CONCLUSIONS 

The major aim of the present study was to elaborate tools, 
which provide comprehensive investigation of remote young 
open star clusters having accurate CCD-photometry. The 
effort was focused on revealing of the Main Sequence turn- 
on areas, important both from the point of young cluster 
dating and from the luminosity and mass function analy- 
sis. To enable the statistical selection of cluster members for 
vast set of open cluster CCD observations with no data on 
adjacent blank fields, and thus, to involve them to luminos- 
ity and mass functions studies we proposed a new approach 
to field star removing technique based on using of galactic 
disk star count model. As we found, due to high confidence 
in interstellar extinction determination, open clusters are 
especially suitable objects for this technique, providing se- 
lection results at least of the same quality level as the stan- 
dard method of the blank field does. The proposed approach 
could be especially valuable for deep observations of young 
clusters, planned to be observed in th e frame of present- 
day surveys (see e.g. lKalirai et al] l200F). which are able not 
only to reach turn-on regions of selected young clusters, but 
also reveal complete H-details down to convective portions 
of luminosity functions. 

The main conclusion from the present study of the 6 
such clusters residing at helio-centric distances of 2-4 kpc is 
that their CCD observations coupled with the Galaxy model 
star counts, wide area statistics provided by all sky cata- 
logues, and theoretical LFs fitted to the observations are 
sufficient to study their population, construction of the de- 
tailed luminosity and mass functions, and age, provided the 
photometry is deep enough to reach the LF H-feature. The 
direct comparison of observed and theoretical LFs instead 
of converting them to stellar mass spectra is of principle im- 
portance in this approach. It provides a standard candle, 
which can be used for a reliable cluster dating. This is very 
difficult otherwise, since the upper MS of young clusters usu- 
ally implemented for this aim, are too steep, and degree of 
stellar evolution is as a rule insignificant to provide reliable 
cluster ages. In fact it just provides an upper estimates for 
the cluster age. The main conclusions can be written as: 

• Contamination of cluster members with field stars is the 
most important factor influencing the lower part of cluster 
CMD and LF. For the clusters under study, field star con- 
tamination varies in the range of 20-50%. 

• The overlapping clusters may considerably distort the 
observed LFs. The "neighbor" contamination, in addition to 
field stars, is equal to 14% in the case of Hogg 15 and to 24% 
in the case of NGC 2383. 

• The LF H-feature was found well beyond statistical 
noise in all the clusters under study. For most of the clus- 
ters its location agrees well with the theoretical prediction. 
We believe, that in the case of NGC 2384, where H-feature 
displays underdeveloped structure, this can be explained by 
complexity of the area, and spatial incompleteness of the 
data. The false H-feature observed in NGC 2383, which cer- 
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tainly is too old to demonstrate it, is due to contamination 
of the cluster field with NGC 2384 stars. 

• Cluster ages derived in the present study are as a rule 
several times lower than those determined from the fitting of 
the theoretical isochrones to the turn off parts of the MS. We 
believe, that our ages are more accurate, than those derived 
from the upper MS. Note that although we used continuous 
star formation model, the derived durations of star forma- 
tion events are only insignificantly differ from zero. 

• Stellar mass spectra of studied clusters are well repre- 
sented with a power law with slopes, which as a rule agree 
with Salpeter' value, within errors. One should keep in mind 
however, that since the observed data concerns the central 
parts of clusters, the derived mass spectra might be different 
from the IMF due to mass segregation effect. The unusually 
flat mass function of NGC 2384 can be treated e.g. as an 
indirect evidence of higher extent of this cluster and of pres- 
ence of mass segregation effect. 
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